Ab initio calculations of H c2 for Nb, NbSe 2 , and MgB 2 
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We report on quantitative calculations of the upper critical field H C 2 for clean type-II supercon- 
ductors Nb, NbSe2, and MgB2 using Fermi surfaces from ab initio electronic structure calculations. 
The results for Nb and NbSe2 excellently reproduce both temperature and directional dependences 
of measured H& curves without any adjustable parameters, including marked upward curvature of 
NbSe2 near T c . As for MgB2, a good fit is obtained for a ir/a gap ratio of ~0.3, which is close to 
the value from a first-principles strong-coupling theory [H. J. Choi et al. Nature, 418 758 (2002)]. 
Our results indicate essential importance of Fermi surface anisotropy for describing H C 2. 



PACS numbers: 74.25.0p,71. 18. +y, 74.25. Jb 

Quantitative descriptions of nature form an integral 
part of physics. In this context, the density-functional 
theory for normal-state electronic structures has substan- 
tially enhanced our understanding on materials and also 
made materials design possible 1]. The purpose of the 
present paper is to extend those calculations in the nor- 
mal state to type-II superconductors in magnetic fields 
to perform systematic calculations of H C 2 . 

The upper critical field H C 2 is one of the most funda- 
mental quantities in type-II superconductors. From the 
early stage of research on type-II superconductors , it 
has been recognized that Fermi surface anisotropy has 
significant effects on H C 2- However, most of the calcula- 
tions performed so far on H C 2 have used simplified model 
Fermi surfaces and/or phenomenological fitting parame- 
ters. Due to this lack of ab-initio-type calculations, our 
understanding on H C 2 remains at a rather primitive level. 
The only exception is the work by Butler on bec Nb Q, 
where he obtained an excellent agreement with experi- 
ments |5( by using Fermi surfaces and electron-phonon 
interactions from his ab initio calculations. Extending 
this calculation to other materials is expected to improve 
our understanding on H C 2 considerably. 

Recently we have derived an H C 2 equation applica- 
ble to low-symmetry crystals, including Fermi surface 
anisotropy, gap anisotropy, impurity scatterings, and 
strong electron-phonon interactions |fj. We here con- 
centrate on the role of Fermi surface anisotropy and 
carry out clean-limit weak-coupling calculations of H C 2 
for Nb, NbSe2, and MgB2. Sauerzopf et al. Q performed 
a careful experiment on H C 2 of high-purity Nb to report 
some discrepancies frornpreceding experiments 0, 0] and 
the theory by Butler H. Also, H r ?. curves of NbSe 2 
HHE3 and MgB 2 Hi E H Hi remains essentially 
unexplained quantitatively. We here treat Fermi surface 
anisotropy without any adjustable parameters by using 
energy band data from first-principles electronic struc- 
ture calculations. 



We have determined H C 2 by requiring that the smallest 
eigenvalue of the following matrix be zero |fj : 
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Here N denotes the Landau level in expanding the pair 
potential and e n = (2n+l)irT is the Matsubara fre- 
quency (/cb = h = I). The function </>(kp) specifies gap 
anisotropy, which is normalized in such a way that the 
Fermi-surface average (4> 2 ) is equal to unity. The matrix 
K = fC(E n , kp, H) is given explicitly in Ref. 0. We need 
to evaluate the average (c/P/C) in Eq. (JJJ appropriately. 

We have obtained Fermi surfaces and velocities from 
electronic structure calculations within the local density 
app roximation (LDA) 0], using the WIEN2k package 
|l5| which is based on the full-potential linear augmented 
plane wave method. The self-consistent calculations have 
been performed by using finite k points. Then we have 
fitted the energy bands by linear combinations of the star 
functions [lfl ITtI lisj to construct Fermi surfaces on finer 
mesh points. Integrations over Fermi surfaces have been 
performed by the linear tetrahedron method 19] . We 
have checked the convergence by increasing k points. 

We shall focus on the temperature and directional de- 
pendences of H C 2 and use the reduced quantities t = T/T c 
and h*(t) = _ dH H 2 C (lff dt \ t _ 1 ■ When it becomes relevant, 
the absolute value of H c2 is fixed by using an experimen- 
tal value at a particular temperature and field direction. 
This procedure is similar in the normal state to incor- 
porating many-body mass enhancement by using a sin- 
gle point of de Haas-van Alphen data. We have fixed 
0(kp) = 1 for Nb and NbSe2, since gap anisotropy may 
not be substantial in these materials. The spin-orbit in- 
teraction have been included self-consistently for NbSe2. 

bec Nb — We first consider temperature dependence 
of H C 2 in Nb. Figure H shows theoretical curves for the 
angle averaged critical field h* (t) in comparison with the 
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FIG. 1: Normalized angle averaged H C 2 for Nb. Open circles 
are the experiment of Ref.^. The solid line is obtained by us- 
ing the Fermi surface from the LDA calculation, whereas the 
dotted line is the Helfand-Werthamer curve for the spherical 
Fermi surface. The dashed line incorporates the difference in 
mass renormalization among different Fermi surface sheets. 



experiment by Sauerzopf et al. (circles) 7] . The solid line 
with h* (0) = 0.96 has been obtained as above by using the 
Fermi surface from the LDA calculation, whereas the dot- 
ted line with ft*(0) =0.73 is the Helfand-Werthamer re- 
sult [2(J for the spherical Fermi surface. Considering that 
there are no adjustable parameters, our result (solid line) 
shows a fairly good agreement with the experiment, im- 
proving on the Helfand-Werthamer curve substantially. 
However, the value h*(0) = 0.96 still falls well below the 
experimental value ft.*(0) = 1.06. 

We attribute this discrepancy to the difference in 
many-body mass enhancement among different Fermi- 
surface sheets. According to Butler's calculation Q, 
strong but k-independent electron-phonon interactions 
cannot cure the discrepancy, since they merely bring an 
overall enhancement of H C 2 but change temperature and 
directional dependences of h* ft) by less than 2%. On the 
other hand, Crabtree et al. |21| found from their analysis 
on de-Haas van Alphen experiments that the electron- 
phonon renormalization factor A„(kp) differs substan- 
tially by more than 30% among the three Fermi-surface 
sheets (n — 1, 2, 3) of Nb, although it does not vary ap- 
preciably within each sheet. Hence we have performed 
another calculation of H C 2 using the scaled Fermi ve- 
locity v„(k F ) = v„(k F )/(l + A„), where v„(k F ) is the 
bare Fermi velocity from the LDA calculation and A„ 
is a sheet-dependent scaling factor taken from Ref. |2lJ. 
This yields the dashed line in Fig. 1 with h*(0) = 1.03, 
showing a considerable improvement on the bare result. 
This fact indicates the importance of band and momen- 
tum dependences of the electron-phonon couplings. It is 



an interesting problem in the future to see whether the 
effects are reproduced naturally through a first-principles 
strong-coupling calculation of H c2 . 

We next discuss angular dependence of H c2 in Nb. 
This dependence for cubic materials has been studied by 
expanding the relative anisotropy H c2 /H C 2 as 0, E2 

H c2 (Cl,t)/H c2 (t) = l+ (2) 

i=4,6,8,--- 

Here fi denotes direction of the field and P;'s are cu- 
bic harmonics |2^| invariant under cubic symmetry op- 
erations and constructed as linear combinations of the 
spherical harmonics. The coefficients ai(t) specify the 
H c2 anisotropy at temperature t. After the work of But- 
ler 0, Sauerzopf et al. Q carried out a very careful and 
detailed experiment on ai(t). Hence it is very interesting 
to see whether ai's can also be reproduced by ab initio 
calculations. We have evaluated a;'s by least square fit- 
tings for H c2 computed over 50 different directions. 

Figure 2 compares calculated ai(t)(l = A, 6, 8, 10) with 
the experiment by Sauerzopf et al. [7j. The agreements 
are good, especially for the dominant a^t). For exam- 
ple, the change in the curvature of a^it) around i«0.6 
is excellently reproduced. This is not the case for the 
second-largest contribution ae(t), however, and the cal- 
culated curve reproduces only 60% of the experimental 
value at t = 0. An improvement is observed by introduc- 
ing the Fermi-surface-dependent renormalization factor 
A„ in the calculation (dashed line), but there still re- 
mains an appreciable deviation. This discrepancy may 
be attributed to the small gap anisotropy which may be 
present in Nb. Indeed, Seidl et al. reported that the 
gap anisotropy has a dominant effect on the term. 
However, their analysis may not be reliable quantita- 
tively, especially near t = 0, since they applied a theory 
by Teichler [24| which incorporate only the lowest-order 
corrections from the gap and Fermi-surface anisotropics. 
The discrepancy may be removed completely by a first- 
principles strong-coupling calculation of H C 2 where the 
gap anisotropy is naturally included. 

NbSe2 — The hexagonal transition metal dichalco- 
genide 2H-NbSe2 has been studied extensively as a pro- 
totype of anisotropic layered superconductors. Nb atoms 
form triangular lattices in the ab plane, which stack along 
the c axis with two Se layers between them. At low tem- 
peratures preceding superconductivity, NbSe2 undergoes 
an incommensurate charge-density-wave (CDW) transi- 
tion. Several groups HUE] studied H c2 of NbSe2 and 
Nbi_ x Ta a; Se2, where H c2 \\ for the magnetic field along 
the ab plane was observed by more than 3 times larger 
than H C 2± for the field along the c axis at T = Q. Another 
characteristic feature is the marked positive curvature of 
H C 2\\ near T c . However, no quantitative calculations us- 
ing detailed Fermi surfaces have been performed yet. 

We have evaluated H C 2 using the Fermi surface from 
the LDA calculation, ignoring any effects of the CDW 
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FIG. 2: Anisotropy coefficients ai of Nb as a function of tem- 
perature for Z = 4, 6, 8, 10. Solid lines are obtained by a bare 
LDA calculation, whereas dashed lines denote results using 
the renormalized Fermi velocity v„(kp) — v„(kp)/(l + A„). 
Open circles are experiments by Sauerzopf et al. Q. 

instability. Calculated energy bands and Fermi surface 
agree well with a previous one |25| . The Fermi surface can 
be classified into two groups. The first one is cylindrical 
sheets, and the other is a small hole pocket around T 
point which looks like a flattened spheroid. 

Figure |31 shows calculated H c2 \\ and H r 2 \ cu rves in 
comparison with a couple of experiments |8|, Il0| . Tem- 
perature dependence of the anisotropy parameter 7 = 
H C 2\\/H C 2± is also shown in the inset. We observe good 
agreements between the theory and experiments for both 
field directions. Especially, the positive curvature of H&h 
has been reproduced excellently. Although this positive 
curvature has been observed commonly in anisotropic 
superconductors [HI, its origin has not been identified 
clearly so far. Thus, we have shown explicitly that Fermi 
surface anisotropy adds quite a variety to H C 2 curves. 
The reduced critical field h*(t) is also very anisotropic, 
reaching h*\(0) = 1.73 and h* ± (0) = 0.68 in our calcula- 
tion. If we ignore the small Fermi surface sheet around 
r point, the agreement becomes worse as shown by the 
dashed line. Thus, the small Fermi surface sheet cannot 
be neglected for the quantitative understanding of H C 2- 

MgB 2 — Superconductivity in MgB2 found by Naga- 
matsu et al. [27j has attracted much attention. Besides its 
high transition temperature ~40K, it has a unique fea- 
ture that the energy gap A differs substantially in mag- 
nitude between a- and 7r-bands, as predicted by theories 
|28l |22j and confirmed by experiments .30, 31, a 
first-principles calculation by Choi et al. |29j has yielded 
Ao.~6.8meV and A,r~1.8meV on the average as T^0. 

Experiments on H C 2 of uniaxial MgB ? single cry stals 
have been reported by several groups [lj, [l^, llM LL-I ■ Mi- 
ranovic et al. [3j| performed a calculation of H C 2 adopting 
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FIG. 3: H C 2\\ and H C 2± of NbSe2 as a function of temper- 
ature. Crosses and open circles are experiments of Ref. @ 
and Ref. [ToL respectively. The solid lines are from the LDA 
calculation. If we ignore the Fermi surface around F point, 
the dashed lines result. The inset plots calculated anisotropy 
ratio H c2 \\/H c2 ±. 

a model Fermi surface with two spheroids, which yielded 
a qualitative agreement with observed H C 2 . We have car- 
ried out a more detailed calculation using the Fermi sur- 
face from the LDA calculation. As seen above, detailed 
Fermi-surface structures are indispensable for the quan- 
titative understanding of H C 2 . As for the gap anisotropy, 
we have incorporated the ratio cv^A^/Ao- at T = as a 
single parameter in the calculation. 

Figure 0] plots calculated H c2 \\ and H C 2± for MgB2 
using three different values of a. When we assume an 
isotropic gap of a = 1, we could not describe the ex- 
perimental anisotropy H c2 \\/H C 2±. The anisotropy of 
0.6 < t < 1.0 has been reproduced most completely for 
the choice a ~ 0.3. This value is close to the one ob- 
tained by Choi et al. [2^| and consistent with the model 
calculation by Miranovic et al. |33j . Thus, the present 
calculation provides an additional support for the exis- 
tence of two energy gaps via a detailed calculation of 
H C 2- The reduced critical field h*(t) is very anisotropic 
as hUO) = 1.64 and h* x (0) = 0.67 for a = 0.3. However, 
calculated H C 2 for a = 0.3 shows larger anisotropy at 
lower temperatures than the experimental one. This dis- 
crepancy may be removed completely by incorporating 
effects neglected here, such as strong electron-phonon in- 
teractions, Pauli paramagnetism, etc. 

In summary, we have evaluated H C 2 of three classic 
type-II superconductors using Fermi surfaces from first- 
principles electronic structure calculations. For Nb and 
NbSe2 whose superconducting gap anisotropy may not be 
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FIG. 4: H C 2 of MgB2 as a function of temperature. Open 
circles are experiments of Ref. Il2l Calculated H C 2 curves for 
A^/A CT = 0.25, 0.3, and 0.35 are shown by dotted, solid and 
dashed lines, respectively. The inset shows a comparison of 
the anisotropy ratio H c2 \\/H C 2±. 



significant, calculated H C 2 curves satisfactory reproduce 
experimental temperature and directional dependences 
of H c2 without any adjustable parameters. As for MgB 2 , 
a good agreement with experiments follows if we choose 
the gap ratio Att/Act^O.S. Thus, the present study has 
clarified unambiguously the necessity of incorporating de- 
tailed Fermi surface structures in the calculation of H C 2 . 
Such calculations will also form a basic starting point to 
discuss other contributions, and we may obtain unique 
information about superconducting gap anisotropy, etc., 
by comparing calculated and experimental H C 2 curves. 
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